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mixing at the detectable level or significant CP violation in the charm system 
may strongly signify the existence of new physics. In view of the large discovery potential 
associated with the fixed target experiments, the P-meson factories and the r-charm 
factories, we make a further study of the phenomenology of D° — D mixing and CP 
violation in neutral D-meson decays. The generic formulas for the time-dependent and 
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Abstract 



time-integrated decay rates of both coherent and incoherent D°D° events are derived, and 
their approximate expressions up to the second order of the mixing parameters x D and 
y D are presented. Explicitly we discuss D° — D° mixing and various CP-violating signals 
in neutral D decays to the semileptonic final states, the hadronic CP eigenstates, the 
hadronic non-CP eigenstates and the CP- forbidden states. A few non-trivial approaches 
to the separate determination of x D and y D and to the demonstration of direct and indirect 
CP asymmetries in the charm sector are suggested. 
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1 Introduction 



It is well known in particle physics that mixing between a neutral P° meson and its CP- 
conjugate counterpart P° can arise, if both of them couple to a subset of virtual and (or) 
real intermediate states. Such mixing effects provide a mechanism whereby interference in 
the transition amplitudes of P° and P° mesons may occur, leading to the possibility of CP 
violation. Determining the magnitude of P° — P° mixing and probing possible CP- violating 
phenomena in the P° — P° system have been a challenging task for particle physicists. To date, 
K° — K° and P° — B% mixing rates have been measured, and the CP-violating signal induced 
by K° — K° mixing has been unambiguously established jjT| . Many sophisticated experimental 
efforts, such as the programs of (p factories, B factories and high-luminosity hadron machines, 
are being made to discover new signals of CP asymmetries beyond the K° — K° system and to 
precisely measure the Kobayashi-Maskawa (KM) matrix elements. 

The study of mixing and CP violation in the Q = +2/3 quark sector, particularly in the D°— 
D° system, is not only complementary to our knowledge on the K° — K° and B° — B° systems, 
but also important for exploring possible new physics that is out of reach of the standard 
model predictions. The rate of D° — D° mixing is commonly measured by two well-defined 
dimensionless parameters, x D and y D , which correspond to the mass and width differences of 
D° and D° mass eigenstates. The latest E691 data of Fermilab fixed target experiments only 
give an upper bound on D° — D° mixing 0: 

r D ~ *d + Vd < 3 . 7xl0 -3 (L1) 

In the standard model the short- distance contribution to D° — D° mixing is via box diagrams 
and its magnitude is expected to be negligibly small (x s ff' ~ 10~ 5 and y% d - < x s ff' 0) . The 
long-distance effect on D° — D° mixing comes mainly from the real intermediate states of SU(3) 
multiplets, such as 

D° «-> 7T7T , ttK , irK , KK «-> D° , (1.2) 

and is possible to be significant if the SU{3) symmetry is badly broken (e.g., x 1 ^ ~ y)^' ~ 
10 -3 — 10 -2 However, the dispersive approach H] and the heavy quark effective theory 
|J seem to favor a much smaller result for the long-distance contribution: x 1 ^' ~ 10 x x^ d ' 
and x X q - ~ x s ^ respectively. Such theoretical discrepancies indicate our poor understanding 
of the dynamics for D° — D° mixing, hence more efforts in both theory and experiments to 
better constrain the mixing rate are desirable. If calculations based on the standard model can 
reliably limit x D and y D to be well below 10~ 2 , then observation of r D at the level of 10~ 4 or so 
will imply the existence of new physics. On the other hand, improved experimental knowledge 
of r D , in particular the relative magnitude of x D and y D , can definitely clarify the ambiguities 
in current theoretical estimates and shed some light on both the dynamics of D° — D° mixing 
and possible sources of new physics beyond the standard model. 
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The phenomenology of CP violation in the D° — D° system was first developed by Bigi and 
Sanda [j/J, and further summarized by Bigi in ref. ||. These works have outlined the main 
features of D° — D° mixing and CP asymmetries anticipated to appear in neutral P-meson 
decays, although many of their formulas and results are approximate or just for the illustrative 
purpose. The theoretical expectations on the magnitudes of various possible effects are also 
sketched in refs. J7|, ||. 

Recent experimental progress, particularly in observing the doubly Cabibbo suppressed 
decay (DCSD) D° — > K + tt~ ||, constraining the D°—D° mixing rate [p], [| and searching for the 
CP asymmetries in D decays to K + K~ etc |L0|, are quite encouraging. Further experimental 



efforts, based mainly on the high- luminosity fixed target facilities ||11|| , the forthcoming P-meson 



factories and the proposed r-charm factories are underway to approach the above physical 
goals. In view of the large discovery potential associated with these experimental programs, a 
further study of the phenomenology of D° — D° mixing and CP violation in the charm system 
is no doubt necessary and important. 

In this paper we shall on the one hand follow the pioneering work of Bigi and Sanda to 
refine upon the phenomenology of D° — D° mixing and CP violation in neutral D decays, 
and on the other hand investigate some specific possibilities to separately determine x D and 
y D as well as to probe various CP- violating signals in the charm sector. A generic formulism 
for the time-dependent and time-integrated decay rates of both coherent and incoherent D°D° 
events are derived, and their approximate expressions up to 0{x 2 D ) and 0{y 2 D ) are presented. 
Systematically but explicitly, we discuss a variety of D° — D° mixing and CP- violating mea- 
surables in neutral D decays to the semileptonic final states, the hadronic CP eigenstates, 
the hadronic non-CP eigenstates and the CP-forbidden states. We show that it is possible to 
determine the relative magnitude of x D and y D through observation of the dilepton events of 
coherent D°D° decays on the ■0(4.14) resonance at a r-charm factory. A model- independent 
constraint on D° — D° mixing can also be obtained by measuring the decay-time distributions 
of D°/D° — > Ks,l + tt° etc. By use of the isospin analysis and current data, we illustrate final- 
state interactions in D — ► KK and their influence on CP violation. The interplay of D° — D° 
mixing and DCSD effects in incoherent D°D° decays to K^ir^ and in coherent D°D° decays 
to both (/ ± X =F , K ± tt t ) and (i^ ± 7r =F , K ± ir T ) states is analyzed in the presence of CP violation 
and final-state interactions. We take a look at two types of CP- forbidden decays at the ^(3.77) 
and -0(4.14) resonances. Finally the possibility to test the AQ = AC rule and CPT symmetry 
in the D° — D° system is briefly discussed. 

This work is organized as follows. In section 2 we derive the generic formulas for coherent and 
incoherent D°D° decays, and then make some analytical approximations for them. Sections 
3, 4, 5 and 6 are devoted to D° — D° mixing and CP violation in neutral D decays to the 
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semileptonic states, the hadronic CP eigenstates, the hadronic non-CP eigenstates and the 
CP-forbidden states respectively, where some distinctive approaches or examples are discussed 
for determining x D and y D or probing possible CP- violating effects. We summarize our main 
results in section 7 with some comments on tests of the AQ = AC rule and CPT symmetry. 



2 Fundamental Formulas 

We first develop a generic formulism for the time-dependent and time-integrated decays of 
neutral D mesons. Considering the smallness of D° — D° mixing indicated by both experimen- 
tal searches and theoretical estimates, we then make some analytical approximations for the 
obtained decay rates up to the accuracy of 0{x 2 D ) and 0{y 2 D ). 

A. Preliminaries 

In the assumption of CPT invariance, the mass eigenstates of D° and D° mesons can be 
written as 



\Dl) 


= P\D°) 


+ q\ 


\Dh) 


= P\D°) 


- q\ 



(2.1) 

in which the subscripts "L" and "H" stand for Light and Heavy respectively, and (p, q) are 
complex mixing parameters. Sometimes it is more convenient to use the notation 



Q 
P 



exp(i20) , (2.2) 



where is a real CP-violating phase in D° — D° mixing. With the help of the conventions 
CP|P°) = \D°) and CP\D°) = \D°), the relations between the CP eigenstates 

m s \jn±m , m s un^m ,, 3) 

and the mass eigenstates \D L ), \D H ) turn out to be 

in \ p + q \n \ + p ~ q \n \ 

I^h> = P -^\D 2 ) + ^\ Dl ) . (2.4) 
The proper-time evolution of an initially (t = 0) pure D° or D° is given as 

\D° phys (t)) = 9 + (t)\D°) + V(t)|P°) , 

\D° phys (t)) = 9 + (t)\D°) + - q 9-(t)\D°) , (2.5) 
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where 



g + (t) = exp 
g_(t) = exp 



lm 



- I i'» f - I *' 



with the definitions 



m T + m H 

m = — ^ ^ 

2 

r L + r H 



cosh 
sinh 

Am 

Ar 



( lAm "f)^ 
( lAm "f)f 



(2.6) 



m H — m 



L J 



r H . 



(2.7) 



Here m L ^ and r L ( H ) are the mass and width of -Dl(h) respectively. Note that the above 
definitions guarantee Am > and Ar > in most cases. Practically, it is more popular to use 
the following two dimensionless parameters for describing D° — D° mixing: 

Am Ar 
2F ' 



x 



D 



Vd 



(2.8) 



Certainly both x D and y D in most cases are positive (or vanishing). 

B. Rates for incoherent D decays 



The transition amplitude of a neutral D meson decaying to a semileptonic or nonleptonic 
state / can be obtained from eq. (2.5) as follows: 



(m\D° phys (t)) = g + (t)A f + ±g_(t)A f , 
(f\n\D° phys (t)) = g + (t)A f + Zg_(t)A f , 



(2.9) 



where Af = (f\H\D°) and Af = (f\H\D°). For convenience, we also define the ratio of these 
two amplitudes: 

Then the time-dependent probabilities of such decay events are expressed as 



R(D° phys (t) ^ f) « \A f \ 2 exp(-Tt)[C y cosh(y D Tt) + C x cos(x D Tt) 

+ S y smh(y D Tt) + S x sm(x D Ft) 

R(D° phys (t) -> /) oc \A f \ 2 ex V (-Tt)[C y cosh(y D Tt) + C x cos(x D Tt) 

+ 5 y smh(y D rt) + S x sm(x D Tt) 



where 



C y = 



1 + IA/l 2 



1- |A f | 2 



Sy = — ReA f , 
S x = -ImA / ; 



(2.10) 



(2.11) 



(2.12) 
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and 



{^Cy i Sy , C x , S x 



(C y , Sy , C x , S X ^j . 



(2.13) 



To obtain the time-independent decay rates, we integrate eq. (2.11) over t G [0, oo) and get 



R(D°^f) oc \A f \ 



-Cy + ; K-C X + TrSy + 



i-vb v i + ^ i-vb v ' i + 4 



1 r _|_ 1 r» _|_ c- 
l - r n l 



y£ " L + x z D ~ ' l-yfj l + x 2 D 
Eqs. (2.11) and (2.14) are the master formulas for incoherent D decays. 



(2.14) 



Following the same way one can calculate the decay rates of -Dp hys and -Dp hys to /, the 
CP-conjugate state of /. To express the relevant formulas in analogy of eqs. (2.11) and (2.14), 
we define Aj = (f\H\D°), Aj = {f\H\D } and 

At 



Pf 



p. 



(2.15) 



Then R(D° phys (t) -> /), R(D° phys (t) -> /) and R(D° phys -> /), R(D° phys -> /) can be written 
out in terms of and Xf. If / is a CP eigenstate (i.e., |/) = CP\f) = ±|/)), then we get 
Aj = ±A f , Af = ±A f , pj = l/pj and Xj = l/Xj. 

C. Rates for coherent D decays 

For a coherent Dp hys Z)p hys pair at rest, its time-dependent wave function can be written as 
[D° phys (K,t)) ® |^ hys (-K,0) + C\D° phys (-K,t)) ® \D° phys (K,t))} , (2.16) 

where K is the three-momentum vector of the D mesons, and C — ± denotes the charge- 
conjugation parity of this coherent system. The formulas for the time evolution of -Dp hys and 
Dp hys mesons have been given in eq. (2.5). Here we consider the case that one of the two D 
mesons (with momentum K) decays to a final state f\ at proper time t\ and the other (with 
— K) to f 2 at t 2 . fi and f 2 may be either hadronic or semileptonic states. The amplitude of 
such a joint decay mode is given by 



A(h,t i; f 2 ,t 2 )c = -^Aj.Aj^clg^g^tJ + Cg^g+ih) 

^h A hCc [g+itjg+fo) + Cg^g^h) 
where A f . = (fi\H\D°) (with % = 1,2), and 

= ^(1 + CX fl X f2 ) , 
(c = ^(A /2 + CX h ) . 



+ 



(2.17) 



(2.18) 
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Here the definition of A/ 2 and A/ 2 is similar to that of A/ in eq. (2.10). After a lengthy 
calculation we obtain the time-dependent decay rate as follows: 

R(fi,h;f 2 ,t 2 )c cc \A fl \ 2 \A h \ 2 exp(-rt+) x 

'(ltd 2 + Kcl 2 ) cosh(y D Tt c ) - 2Re (C c Cc) smh(y D Tt 



a 



where 



has been defined. 



- (I£d 2 " Kcl 2 ) cos(x D H c ) + 2Im (Cc(c) sm(x D Tt c ) 

t c = h + Ch 



(2.19) 
(2.20) 



The time-independent decay rate is obtainable from eq. (2.19) after the integration of 
R(fi,h; f 2 , h)c over t x E [0, oo) and t 2 E [0, oo): 

1 + (\t- l2 , „ l2 \ 2(l + C)y D . 



R{hJ 2 )o oc I^JVj 



1 C x 



i - vl) 2 

(ied 2 -ic 



(led 2 + led 



(i - z/i)) 2 



'Re (^cCc) 



+ 



2{l + C)x D 
(1 + x 2 ^ 2 



Im (&Cc) 



(2.21) 



(1 + X D ) ■ V x -r ^ D ) 

We see that two interference terms Re(£££c) and Im(^Cc) disappear in the case of C = — 1, 
independent of the final states /i and f 2 . 

In a similar way, one can calculate the joint decay rates of (-Dphys-^phys)c to (fif 2 ), (fif 2 ) 
or (/1/2), where /1 and /2 are CP-conjugate states of fi and f 2 respectively. 

D. Analytical approximations 

In the standard model, the magnitudes of x D and y D are expected to be very small, at most 
of the order 10 -2 (see, e.g., refs. [[§, [5], ||). The current experimental constraints on D° — D° 
mixing give x 2 D + y 2 D < 7.4 x 10~ 3 (see eq. (1.1)), which implies x D < 0.086 and y D < 0.086. 
Due to the smallness of x D and y D , the generic formulas obtained above can be approximately 
simplified to a good degree of accuracy. 

Up to the accuracy of 0(x 2 D ) and 0{y 2 D ) for every distinctive term, the time-dependent 
decay rates in eq. (2.11) are approximated as 

1 

4 



R(D° phyB (t) ^ f) oc \A f \ 2 exp(-Tt) 1 + j(x 2 D + yl) | X f \ 2 T 2 t 

1 



R(D° phys (t) - /) cx \Aj 



P 

q 

- \ x 



- \ [ x d ~ Vd) ^ ~ (VD ReX f + x D lm\ f )Tt 
2 



exp(— Ft) 



' D 



|A/| 2 + ~(x 2 D + yl)rH 2 



y 2 D )\X f \ 2 T 2 t 2 - (y D Re\ f - x D lm\ f )Tt 



(2.22) 
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Similarly we obtain the approximate decay rates for D® hys (t) — > / and D® hys (t) — ► /: 

2 



^PphysW 



14 -I 2 


g 


P 






-si 


x 2 D 



exp(— Ft) 



|A 



f\ 



+ - d & + vl) r 2 t 2 



vl) IAHW 



4 



R(D° ph Jt) ^ f) cc |A / fexp(-R) 



4 



x 2 D —y 2 D )v 2 t 2 — (y D Re\f + x^ImAj) 



(2.23) 



The time-independent rates for these four processes turn out to be 

7 A' \ -s. I I I - 1 ' 

phys 



R(D° 



f) oc |^| 



/I 



i?(.D 



phys 



/) oc \Aj 



i + o ( x2 D+y 2 D)\hf 



o ( z D - !/£> ) - ( 2/z? ReA /~ + a^ImA; 



(2.24) 



and 



R(D° 



phys 



/) oc \Aj 



\\ f \ 2 



- i; ( x d ~ Vd) - (VD ReX f ~ x D lmX f ) 



R(D° 



phys 



/) oc \A f \ 



Q 
p 



am 



y D Re\f — x D lm\f 



(2.25) 



The formulas listed above are very useful for the study of neutral D decays in fixed target 
experiments or at .B-meson factories. Here no assumption has been made for the magnitudes 
of | A/ 1 and \Xf\- If they are considerably smaller than unity, e.g., in the DCSDs, then much 
simpler expressions can be drawn from eqs. (2.22) to (2.25). 

It is common knowledge that the decay-time distributions of coherent (-Dphys-Dphys)c* pairs 



cannot be measured at a symmetric e + e~ collider [14]. Since the presently-proposed r-charm 
factories are all based on symmetric e + e~ colisions, it is more practical to study the time- 
integrated decays of (D® hys D® hys )c pairs. For completeness we shall present some important 
formulas for the decay-time distributions of (D^ hys Dp hys ) c events, with the assumption of an 
asymmetric r-charm factory, in Appendix A. Such a work might be of purely academic sense, 
but it could also be useful in the future experiments of charm physics. 

In the approximations up to 0(x 2 D ) and 0{y 2 D ), the time-integrated rates for (-Dphys-^phys)c 
decaying coherently to (/1/2), (/1/2), (/1/2) and (/1/2) states are obtained from eq. (2.21) as 
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follows: 



and 



R(hJ 2 )c oc \A fl \ 2 \A h \ 



+ 



R{hJ 2 )c oc \A fl \ 2 \A f , 



{(2 + C) (x 2 D + y 2 D )\l + C\ fl \ h \ 2 

2-(2 + C) (x 2 D -y 2 D )] \\ h +C\ h \ 2 
- 2(1 + C)y D [(l + lA^I 2 ) ReA /2 + C (l + |A /2 | 2 ) ReA^j 
+ 2(1 + C)x D [(l - |A /X | 2 ) ImA /2 + C (l - |A /2 | 2 ) ImA/J } 

{(2 + C) (4 + 2/ 2 ,)|l + CA /l A /2 | 2 



+ 



2-(2 + C) 4 



|A /2 + CA 7l 



-2(l + C)y D 
+ 2(l + C)x D 



l + |A /l | 2 )ReA /2 + C(l + |A /2 | 2 )ReA j 
(l-|A /l | 2 )lmA /2 + C(l-|A /2 | 2 )lmA /l ]} ; (2.26) 



R(fuh)c oc |A /l | 2 |%| 2 {(2 + C)(4 + y 2 ) )|A /2 +CA /l | 



+ 



2-( 2 + C)(xj D -yf ) \l + CX fl X f , 
2(1 + C)y D 
2(l + C)x D 



l + |A A | 2 )lteA J 5 l + C7(l + |A A | 2 ) ReA A ] 
1-|A/J 2 ) ImA /2 +C(l-|A /2 | 2 ) ImAj} 



R(fuf2)c oc |%| 2 |A /2 | 2 {(2 + C)(x 2 3 + ^)|A /2 +CA /l | 2 
+ [2-(2 + C)(x 2 D -y 2 D )}\l + CX fl X h \ 2 

- 2(1 + C)y D [(l + lA^I 2 ) ReA /2 + C (l + |A /2 | 2 ) ReA 7l 

- 2(1 + C)x D [(l - |A 7l | 2 ) ImA /2 + C (l - |A /2 | 2 ) ImA 7l ] } . (2.27) 

Taking /i = K + l^u l or /i = K~l Jr v l for example, eqs. (2.26) and (2.27) can be simplified sig- 
nificantly. Such semileptonic decay modes, which are flavor-specific, play the role in identifying 
the flavor of the other D meson decaying to f 2 or f 2 . 



3 Semileptonic D decays 

The manifestation of D° — D° mixing and CP violation in the semileptonic decays of neutral 
D mesons is relatively simple, since such transitions are flavor-specific in the standard model 
or some of its extensions. Due to the flavor specification of D° — > l + X~ and D° — > l~X + , it is 
not necessary to study the time dependence of -Dp hys and -Dp hys decay modes. 

A. D° - D° mixing and CP violation 
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For fixed target experiments or e + e~ collisions at the T(4S) resonance, the produced D° 
and D° mesons are incoherent. Knowledge of D° — D° mixing is expected to come from ratios 
of the wrong-sign to right-sign events of semileptonic D decays: 

R(D° phys - l + X-) ' R(D° phys - 1-X+) ■ 1 • 1 

By use of eq. (2.14), we find 



2 1 
1 — a 

1 + a 



■ 



a 



1 + a 



(3.2) 



where a = (1 — y 2 D )/{l + x 2 D ). Note that \q/p\ ^ 1 signifies CP violation in D° — D° mixing. 
To fit more accurate data in the near future, we prefer the following mixing parameter: 

_ r + f 1 - a 

r ° = ~Y~ = w TT^ (3 - 3) 

with w = (\q/p\ 2 + \p/q\ 2 )/2. For \q/p\ — 1 ~ ±1%, the value of w deviates less than 0.1% from 
unity. Thus this overall factor of r D is safely negligible. In the approximation of x D <C 1 and 
j/ D < 1, one obtains 

r D - ■ (3.4) 

The latest E691 data [§] give r ps f m r D < 0.37% for small x D and y D , where \q/p\ ~ \p/q\ ~ 1, 
a worse approximation than w « 1, has been used. 

The CP asymmetry between a semileptonic decay mode and its CP-conjugate counterpart 
is defined as 



A D = 



R(D° phys -» 1+X-) - R(D° phys -» /-X+) 



Straightforwardly, we get 



A ^ = mLm! > A D = . (3.6) 



If Ad is at the level of 10 3 or so, it can be measured to three standard deviations with about 
10 7 wrong-sign events. 

It should be noted that the asymmetry A & may be nonvanishing if there exists new physics 
affecting the semileptonic D decays. For example, either the violation of CPT symmetry or 
that of the AQ = AC rule can lead to A^ ^ 0. Even if the AQ = AC rule and CPT invariance 
hold, A D ^ is still possible in consequence of the phase shifts from final-state electromagnetic 
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interactions or the CP- violating contributions of non-standard electroweak models to the tree- 
level processes under discussion. Hence all such fine effects should be kept in mind and carefully 
evaluated when one wants to isolate one of them from the others. 



As pointed out by Bigi in ref. ||, a nonvanishing value for r D might only be a secondary 
signature of D° — D° mixing, because the presence of AQ = —AC transitions would contribute 
to r D in a significant and time-independent way. For the purpose of illustration, we shall 
specifically calculate this effect on the magnitudes of r D and in the following. 

B. Effect of AQ = —AC transitions on r D and Ad 

Within the standard model the processes D° — > l~X + and D° — > l + X~ are forbidden 
according to the AQ = AC rule. New physics beyond the standard model may allow AQ = 
—AC transitions, which affect the parameters of D° — D° mixing and CP violation. In the 
assumption of CPT symmetry and the neglect of final-state electromagnetic interactions, the 
decay amplitudes of D° and D° to Z ± X =F can be factorized as follows: 



(i + x-\n\D°) 
(rx + \H\D°) 



A t , (l+X-\H\D°) = a t Ar, 
A* , (l-X + \H\D ) = a*A* , 



(3.7) 



where a l measures the AQ = —AC transition amplitude. With the help of eq. (2.14) and 
notations 



we obtain 



\ - q 
V 



V 

-a 

q 



i ■ 



(3.8) 



R(D 
R(D° 



o 

phys 



phys 



l + X~) oc |^| 2 
rx + ) oc lAA 2 



'1 + a) + (1 — ct) | A + 1 2 — 2y D Re\ + — 2ax D \m.\ A 
l + a) + (l-a)|A_| 2 -2?/ 1) ReA_-2ax D ImA_ ; (3.9) 



and 



R(D 



o 

phys 



R(D° 



phys 



l + X~) oc \A\ 2 
rX + ) oc \A L \ 2 



1 -a) + (l + a) I A. 



2y D Re\ + + 2ax D lm\. 



[1 - a) + (1 + a)|A_| 2 - 2y D ReA_ + 2ax D ImA_ . (3.10) 



For small \a l \ (e.g., \a l \ ~ x D or y D ), the original mixing parameters r and r take the following 
forms: 



r w r + \cr, 



■2y 



D 



1 + a 



ReA. 



2ax 



f > « f+\a,\ 2 -^-ReX 



1 + a 

2ax 



^ImA 



+ 5 



1 + a 



1 + a 



^ImA_ . 



(3.11) 
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As a consequence, 



I 

r D = 



In two extreme cases a. 



= and r 



y_p 

l + a 



Re(A+ + A_) 



— ^Im(A+ + A_). (3.12) 
l + a 



D 



0, we obtain r', 



D 



r D and r' D 



a l | 2 respectively. This 



implies that a nonzero value for r' D might not result exclusively from D° — D° mixing. For this 
reason, the study of D° — D° mixing in some other decay modes of neutral D mesons (e.g., 
D°/D° — > K ± tc t ) is necessary in order to pin down possible new physics in the charm sector. 

The magnitudes of CP asymmetries Ad and A^ might be affected by the AQ = —AC 
transitions too. In the approximation of \cx l | <C 1 and r D < 1, we find that becomes 

ax D 



A' 



D 



l A 



D 



Vd 



■Re(A_ - A+) 



-Im(A_ - A+) . 



(3.13) 



l+a l+a 
Note that nonvanishing A^, comes from the interference between the D° — D° mixing and 
AQ = —AC amplitudes, i.e., either r D = or <r l =0 can give rise to A' D = A D = 0. If A^ = 
is assumed, then one obtains A' D ps x D ImA + . Since both x D and \a l \ are expected to be very 
small (even vanishing), observation of the CP asymmetry A' D may be practically impossible. 

C. Separate determination of x D and y D 

Current theoretical estimates for the sizes of x D and y D have dramatical discrepancies due 
to the difficulty in dealing with the long-distance interactions 0, [!| ||. Hence a separate 
determination of these two mixing parameters from direct measurements is very necessary 
(T5|, |ll|] . Here we propose a time-independent method to probe the relative size of x D and y D in 
the dilepton events of coherent Dp hys /)p hys decays at the -0(4.14) resonance. In our calculations 
both CPT invariance and the AQ = AC rule are assumed to hold exactly. 

For a r-charm factory running at the ^(4.14) resonance, the coherent D°D° events can be 
produced through ^(4.14) — >■ 'j(D°D ) c =+ or ^(4.14) -» tt (D°D ) c =-, where C stands for the 



charge-conjugation parity [|nj. The generic formulas for the joint decay rates of two D mesons 
have been given in eq. (2.21). For our present purpose, we only consider the primary dilepton 
events which are directly emitted from the coherent (-Dphys-^phys)c decays. Let and 
denote the time-integrated numbers of like-sign and opposite-sign dilepton events, respectively. 
By use of eq. (2.21), we obtain 



c 



N, 



c 



N, 



c 



N, 



c 



p 


2 - 


Q 






2 r 


Q 




P 






N, 



c 



2N, 



c 



(i - Vl f 



vl) 2 



+ 



1 C X 

2T2 



X 



D) 



(3.14) 
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where Nc is the normalization factor proportional to the rates of semileptonic D° and D° 
decays. It is easy to check that the relation 

N+ + N~ = N+ + NI- (3.15) 

holds stringently, and it is independent of the magnitudes of D° — D° mixing and CP violation. 

Of course a coherent D°D° pair with C = — can be straightforwardly produced from the 
decay of the r/>(3.77) resonance. Its time-independent decay rates to the like-sign and opposite- 
sign dileptons obey eq. (3.14) too. At a r-charm factory the (D°D°)c=- decays at both the 
^(3.77) and -0(4.14) resonances will be measured, and a combination of them might increase 
the sensitiveness of our approach to probing D° — D° mixing. 

Usually one is interested in the following two types of observables: 

a r = -^rr — , r r = c , c , 3.16 

which signify nonvanishing CP violation and D° — D° mixing, respectively. Explicitly, we find 

a - ~ a + ~ ^ D ~ ■ (3.17) 

If a_ or a + is of the order 10 -3 , it can be measured to three standard deviations at the second- 
round experiments of a r-charm factory with about 10 7 like-sign dileptons (or equivalently, 
about 10 10 D°D° events). Furthermore, 

1-a (3 -a 2 

r_ = w , r , = w , (3.18) 

1 + a + 13 + a 2 v ' 

where j3 = (1 + y 2 D )/{l —x 2 D ). One can see that r_ = r D holds without any approximation. For 
small x D and y D , we have 

2 2 

These two approximate results have well been known in the literature (see, e.g., refs. [7j U). 
In such an approximation, however, the relative size of x 2 D and y 2 D cannot be determined. 

To distinguish between the different contributions of x D and y D to D° — D° mixing, one 
has to measure r ± as precisely as possible. With the help of eq. (3.18), we show that the 
magnitudes of x D and y D can be separately determined as follows: 

2 fl + r_ l + 3r_ l + r + \/l+r_ l+r + \ _1 

D \l — r_ l — r_ \ — r + )\\ — r_ l — r + J 

Vl = fl^.i^ _ hzl±) (hzl± - i^V 1 . (3.20) 



13 



Here it is worth emphasizing that w as the overall (and common) factor of r D , r_ and r + can 
be safely neglected. In the approximations up to 0(r^_) and 0(r+), we obtain two simpler 
relations: 

xl-yl^2 T -±^, xl + yl « Ar_ r -±^. (3.21) 

r + — r_ r + — r_ 

Thus it is crucial to examine the deviation of the ratio r + /r_ from 3, in order to find the 
difference between x 2 D and y 2 D . Instructively, we consider three special cases [pi: 



T 

x d » Vd =^ — « 3 + 2r_ > 3 

r_ 



-d ~ Vd ~ 3 - 9rl « 3 , 



x D « y D ^ ^ 3-2r_ < 3. (3.22) 

These relations can be directly derived from eq. (3.18) or (3.20). If r_ is close to the current 
experimental bound (i.e., r _ = r D ps (x|, + y 2 D )/2 < 0.37%), then measurements of r + /r_ to 
the accuracy of 10~ 4 can definitely establish the relative magnitude of x D and y D . To this goal, 



about 10 8 like-sign dileptons (or equivalently, about 10 11 events of (D°D°)c=- and (D°D°) 
pairs) are needed. 

For illustration, we take a look at the changes of the measurable 

7 = — - 3 (3.23) 

with x D by fixing the value of y D . Allowing 10~ 4 < r < 3.7 x 10~ 3 and taking y D = 0.001, 
0.04 and 0.08 respectively, we plot 7 as the function of x D in fig. 1. It is clear that 7 reflects 
the information about the relative magnitude of x D and y D , and it can be detected if r_ is of 
the order 10~ 3 or so. 

In the assumption of a dedicated accelerator running for one year at an average luminosity 
of 10 33 s _1 cm~ 2 , about 10 7 events of / y(D°D°)c=+ and the similar number of 7T (D D°)c=- 



are expected to be produced at the -0(4.14) resonance [12]]. The precision of 10 to 10 in 
measurements of r_ and r + is achievable if one assumes zero background and enough running 
time Hl2| , |i~7H , and then the similar precision can be obtained for the ratio r + /r_ without much 
more experimental effort (see eq. (3.22) for illustration). If D° — D° mixing were at the 
level of r D ~ 10 -3 (or at least r D > 10~ 4 ), then the relative magnitude of x D and y D should 
be detectable in the second-round experiments of a r-charm factory (beyond the one under 
consideration at present). 
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4 Neutral D decays to CP eigenstates 



Neutral D-meson decays to hadronic CP eigenstates / (i.e., |/) = CP\f) = ±|/)), such as 
/ = 7T + 7r _ and KsiT , are of particular interest for the study of CP violation in the charm 
sector. The formulas for their decay rates derived in section 2 can be simplified because of the 
relations Aj = ±Af, Aj = ±Af, pj = l/pf and Xf = 1/A/. If one takes \q/p\ = 1 in some 
cases, then A j = \*j is obtainable ]18| . 



A. Three sources of CP violation 

In the experimental analyses of incoherent D decays, the combined time- dependent rates 

K ± (t) = R(D° phys (t) - /) ± R(D° phys (t) - /) (4.1) 

are commonly used. For convenience in expressing our analytical results, we first define 

Vr „ _ "21mA; MJ _ 2ReA/ 



U 



1-IA,' 2 



1 + IA/l 2 



V 



/ 



1 + IA/l 2 



/ 



1 + IA/l 2 



which satisfy a concise sum rule 

U ) + V f + W f = 1 • 
With the help of eq. (2.11), we obtain 

Tl + (t) = 72oexp(-n) \cosh(y D Vt) - W f smh(y D Vt) 
—ApUf cos(x D Tt) — ApVf sxn(x D Tt) , 

H-(t) = Kq exp(-H) \-A D cosh(y D Tt) + A D W f smk{y D Tt) 
+Uf cos{x D Tt) + Vj sin(x D rt) , 

where 

oc ^l^/| 2 (l + |^| 2 ) (l + IA/l 2 ) 

is a normalization factor, and Ad = (\p\ 2 — |g| 2 )/(|p| 2 + \q\ 2 ) is related to Ad through 

2A D 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



A D 



i + A 2 



(4.6) 



D 



To properly describe the signal of direct CP violation in neutral D decays, we further define 

1 ~ IP/I 2 



T/ " 1 + IP/I 2 ' 
By use of eq. (2.10), we obtain the relation between Tf and Uf. 

T f + A D 



(4.7) 



U 



f 



1 + A D T f 
15 



(4.8) 



It is clear that Ad, 7} and V/ measure the CP asymmetry in D° — D° mixing, the direct CP 
asymmetry in the transition amplitudes of D decays, and the indirect CP asymmetry arising 
from the interplay of decay and D° — D° mixing, respectively. These sources of CP- violating 
effects appear in TZ±(t) simultaneously, but they have different time distributions and can in 
principle be distinguished from one another |jl9| . The magnitudes of A^, 7} and V/ are expected 



to be very small (e.g., at the percent level in some extensions of the standard electroweak model 
|20||). In contrast, the CP-conserving quantity W/ should be of 0(1). Thus the cos(x D Tt) and 
sin(x D rt) terms are considerably suppressed in 1Z + {t). This interesting feature implies that the 
mixing parameter y D is possible to be constrained from the measurement of the flavor-untagged 
decay rate TZ + (t). We shall discuss this possibility for some neutral D-meson decays in the next 
subsection. 

In lowest-order approximations, we keep only the leading terms of A^, 7} and V/ in 7Z±(t). 
Then the CP- violating observable is given as 

At) = « -Ad + U, + x D V f Vt ' x D V f Vt . (4.9) 

One can see that has little contribution to A(t), and the term Tf is almost independent of 
the decay time t. 

Integrating 7Z±(t) over t G [0, oo), we obtain the time-independent decay rates as follows: 

n 



i-vb 
n 



1 - y D W f - aA D (U f + x D V f ) 

(y D Wf-l)A D + a(U f + x D V f ) , (4.10) 



where a has been given in section 3- A. The corresponding CP asymmetry turns out to be 

A = |- « -A D + U f + x D V f w T f + x D V f (4.11) 
in the leading-order approximation. 

At the t/>(3.77) and ip(4.l4) resonances, the produced D°D° pair may exist in a coherent 
state until one of them decays. Hence we can use the semileptonic decay of one D meson to 
tag the flavor of the other meson decaying to a flavor-nonspecific CP eigenstate /. The time- 
integrated rates of such joint decays can be read off from eq. (2.21). We are more interested in 
the following combinations of decay rates: 

Ot(C) ee R(l-X + ,f) c ± R{l+X-J) c . (4.12) 

After some straightforward calculations, we obtain 

Q + (C) = Q \(l +Cy 2 D ) - (l - Cx 2 D ) « 2 Ad U f 
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- (1 + C) (y D W f + x D a 2 A D V f 
fi_(C) = n [-(l + Cy 2 D )A D + (l-Cx 2 D )a 2 U f 



A- = 7 f\ 7 , \ ~ -Ad + U f « 7} 



+ (1 + CO (i/oA^W/ + ^a 2 V/)J , (4.13) 

where 

fi oc 2 ^° |A|2 2 , (4.14) 

and other quantities have been defined before. Keeping the leading terms of A D , 7} and V/, 
we get the CP asymmetries for C = ± follows: 

ft_(C 
Q + (C 

A+ = n+(c = +) w " Ad + W/ + 2XdV/ w T/ + 2x ° Vf ■ (4,15) 

Indeed, A- is exactly independent of the indirect CP- violating term V/. The asymmetry A+ 
is mainly composed of two sources of CP violation. Comparing eq. (4.15) with eq. (4.11), one 
can see that there exists an interesting relation among three time- independent CP measurables: 

A- + A + w 2A. (4.16) 

This result should be testable in a variety of neutral D decays to CP eigenstates. 

B. An approach to constrain y D and x D 

It has been pointed out that y D might be probed through measurements of the singly 
Cabibbo-suppressed decays D® hys (t) — > K + K~ and 7r + 7r~, if CP conservation could hold in 
them [|Tj|. This idea can be straightforwardly understood from the combined decay rates 7Z+(t) 
in eq. (4.4). Assuming CP invariance, i.e., A j d = 7j = V/ = and W/ = 1 (or W/ = —1 for 
CP-odd final states), we find 

K+(t) = TZ Q exp(-rt) [cosh^rt) - smh(y D Tt)\ = 7Z exp [-(1 + y D )Tt] (4.17) 

with TZq oc 2|A/| 2 . Because (1 + y D )T = Tl, the signature of D° — D° mixing is indeed 
a deviation of the slope of lZ+(t) from exp(— Ft). Since T can be measured via other ap- 
proaches, one is then able to constrain the magnitude of y D . The above method depends 
strongly upon the assumption of CP conservation in D decays, hence it may not work well in 
practice. Subsequently we shall show that a model-independent constraint on y D (or x D ) is in- 
deed achievable, without any special assumption, from measuring the decay-time distributions 
ofD° phy s(t)/D° phys (t)^K SiL + n etc. 

Exactly speaking, Ks^° and Kltt are not CP eigenstates due to the existence of small CP 
violation in K° — K° mixing. Here we want to keep this CP-violating contribution to D decays 
(measured by e), but it can be safely neglected in most cases. 
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In the standard model the transitions D° — ► K°7r° and D° — ► i^°7r° (and their CP-conjugate 
processes) are Cabibbo-allowed and doubly Cabibbo-suppressed, respectively. Both of them 
occur only through the tree-level quark diagrams, as illustrated in fig. 2. Since any new physics 
cannot significantly affect the direct decays of charm quark via the tree-level VT-mediated 
graphs PH| , one expects that fig. 2 remains to be a valid quark-diagram description of the 
above-mentioned decay modes even beyond the standard model. Indeed significant new physics 
may exist in D° — D° mixing and the loop-induced penguin transitions of D mesons . 
The processes D° — > K$,l + ^ an d D° — > K$,l + 71-0 take place through fig. 2 with K° — K° 
mixing in the final states. The mass eigenstates of K° and K° mesons can be written as f\ 

\K 8 ) = (l + e)\K°) + (l-e)\K°), 

\K L ) = (l + e)\K°) - (l-e)\K°), (4.18) 

where the complex parameter e has been unambiguously measured (|e| w 2.27 x 10 -3 and 
4> t ~ 43.6° IQ). Note again that we do not assume Ks^° and K^tc to be the exact CP 
eigenstates, although such an assumption is safely allowed by our main results presented later 
on. The overall decay amplitudes of D / D° — > Ks,l + tt° are then given by 

A{D^K s ,l + A = (l + e*)A KOw0 ± (1-e*)^, 

A(D° -> K S;L + tt°) = (l + e*)A w ± (1 - e*) A ROn0 . (4.19) 

Here Ak°k° etc can be factorized as follows: 

A K07T o = (VcdV^Tx exp(Wi) , Ago^o = (V cs V* d )T 2 exp(i5 2 ) ; 

Arw = exp(Wi) , A K o T o = (V*V ud )T 2 exp(iS 2 ) , (4.20) 

where Kjs etc are the KM matrix elements, T% and T 2 stand for the real (positive) hadronic 
matrix elements, S\ and 5 2 are the corresponding strong phases. Denoting h = T 2 jT\ and 
<5 = <$2 — 5i, we obtain 

:i + ^)(K;K rf )/^exp(W) + (1 - e*)(V; d V us ) 



PKstt — + 



(l-e*)(V cs V u * d )hexp(i5) + (I + e *)(V cd V* 



(l + e*)(V; s V ud )hexp(i5) - (1 - e *)(V; d V us ) 
PKl7I ° " (l-£*W cs V: d )hexp(i5) - (l + e*)(V cd V u * s )- { ' ) 



By use of the Wolfenstein parameter A « 0.22, we have |e| « A 4 , V* s V u d « 1 and V^Ks ~ —A 



2 



Furthermore, /i ~ 1 is anticipated in the factorization approximation f24j . As a consequence, 

P^ s .o ~ 1 + 2e* , Pkl7t0 « -p w (4.22) 

hold to a good degree of accuracy. This result implies that the direct CP asymmetries in 
D° hys (t) I D° hya (t) -> K 5jL + vr° are dominated by K° - K° mixing [0: 

T Ks7r o « 7^ £7r o « -2Ree « -2|e|cos0 e . (4.23) 



2 For simplicity, we neglect the common normalization factor l/y2(l + |e| 2 ) for |jfs) and 
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Explicitly, we get T Ksn o m T Kl7T o ps —3.3 x 10 3 . 

For simplicity, we shall use the notation q/p = \q/p\ exp(i20) (see eq. (2.2)) later on. With 
the help of eqs. (4.22) and (4.23) as well as the reasonable assumption \A D \ < 10~ 2 , we obtain 





w +U Kl7T o 


=a Ad — 2 e cos0 e , 




V KS 7T° 


~ - v Kl ^ 


=» 2 e sin e cos(20) — 


sin(20) , 


Wk s *° 




w 2 e sin e sin (20) + 


cos(20) 



(4.24) 

in good approximations. Clearly the unknown new physics may enter Vk s tt° an d Wkstt through 
the D° — D° mixing phase 0. Within the standard model one expects ~ 0, thus Vk s tt° ~ 
3.1 x 10~ 3 and Wk s tt° ~ 1- Beyond the standard model it is possible that the magnitudes of 
Vk s tt° an d y^Ks-K are dominated by sin(20) and cos(20), respectively. The quantities Vk l -k° 
and Wk l -k° are m the similar situation. 

Due to the smallness of x D and y D , some analytical approximations can be made for lZ±(t) 
in. eq. (4.4) up to 0(x D ) and 0(y D ). Taking eqs. (4.23) and (4.24) into account, we find 

n^\t) ~ k$ s *° exp(-n) [i + xr 2 t 2 - m] , 

n+ Ln °(t) ~ 7^ L "° exp(-rt) [1 + xr 2 t 2 + m] , (4.25) 

where X and F are functions of x D and y D : 



X ~ — 
2 



y|, + z|>A D (A d - 2 1 e I cos 



F ~ 21 el sin 1 



'e 



y D sin(20) + iflAfl cos(20) + y D cos(20) — x d Ad sin(20) . (4.26) 



We can see that X and F vanish in the absence of D° — D° mixing, and the contribution of x D 
to them is significantly suppressed by Ad- Naively one might expect to measure the deviations 
of lZ+ sn (t) and 7Z+ L7T (t) from 7Zq s7T exp(— Ft) and 7Zq l7T exp(— Ft), respectively, in order to 
determine the sizes of X and F. However this is very difficult, if not even practically impossible, 
because of the smallness of X and F. The interesting point here is that a comparison between 
the time distributions of 7l+ s7r (t) and K+ L7T (t) can definitely constrain the magnitude of F. 
In view of |e| ~ 1CT 3 , x D < 0.086, y D < 0.086 and \A D \ < 10~ 2 from our present experimental 
knowledge (and theoretical expectation), only the y D cos (20) term of F is possible to be at the 
percent level (magnitudes of the other three terms in F are all below 10~ 3 ). One can conclude 
that the detectable signal of F has to be at the percent level and it must come mainly from 
the width difference of D° and D° mass eigenstates. For illustration, the time distributions of 
TZ+ sn °(t) and n+ L7T ° (t) are depicted in fig. 3 by taking y D ps 0.08 and ps 0. We see that 
around Ft — 2 the difference between 1Z+ S7T (t)/lZ^ s7T and 1Z^ L1T (t)/7l^ L7T can be as large as 
5%, allowing us to extract a signal of D° — D° mixing provided that the accuracy of practical 
measurements is good enough. 
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The asymmetry between 7Z+ sn (£) and TZ^. L1T (t) can be given as 



A LS {t) 



2A 2 ^(i + xr¥ 



YTt . 



(4.27) 



nl 8 *\t) + n% L * \t) " h 

Indeed the coefficient — 2A 2 cosS/h measures the decay-rate asymmetry between D° — > Ks + tc° 
or their flavor-conjugate processes) |19], j24|. The measurement of Als{P) 



and -D c 



Kr 



7T 



allows us to extract the magnitude of K. To give one a numerical feeling, the changes of Aisit) 
with t are illustrated in fig. 4 by assuming X w yjj/2, Y m y D cos (20), ft, w 1, 5 ~ and taking 
y D ~ 0.08, | cos(20)| « 1. It is clear that a large signal of y D should be detectable from Als^)- 

The effects of D° — D° mixing and CP violation also manifest themselves in the combined 
rates Tl^\t) and TZ Kl7T " (t): 



1Z Ks7T °(t) « 7^ s7r ° exp(-rt) 
U K ^\t) « ^^°exp(-n) 



-2|e| cos < 
-2|e| cos i 



+ XT 2 t 2 + YTt 
+ XT 2 t 2 - Y'Tt 



(4.28) 



where 



x D \e\ cos ( 
2 1 el sin 6 e 



y D Az)sin(20) +x D cos(20)| + y D A D cos(2<f)) - x D sin(20) (4.29) 



with r D « (x 2 D + y1))/2. Obviously the CP asymmetry induced by K° — K° mixing (i.e., Ree) 
plays an important role in the decay modes under discussion. The contribution of CP violation 
in D° — D° mixing (i.e., A^) to Y' is not significant even if <p ~ 0. If new physics considerably 
enhances x D and 0, e.g., x D ~ 10~ 2 and |sin(20)| ~ 1, then 1Z K _ S ' K (t) and 7l^ L7T (t) will be 
dominated by the CP asymmetry arising from the interplay of decay and D° — D° mixing. In 
other words, the signals of CP asymmetries 

n« s7T °(t) 



AK S TT<i(t) 
A Kl 7T° (t) 



K 
K 



K s n° 
+ 

K L n° 



(t) 



(t) 



-2 e cos< 



-2|e| cosi 



XT 2 t 2 + Y'Tt 



+ x'tH 2 



Y'Tt 



(4.30) 



at the percent level will indicate that new physics is definitely present in D° — D° mixing (e.g., 
| sin(20)| > 0.5) and the magnitude of x D must be of order 10 2 . Taking x D « 0.08, « tt/4 
and Ap « for example, we illustrate the time distributions of 7^5 s7r (t) and 7Z^: L7T (t) in fig. 
5. We find that around Tt = 1 the magnitudes of the decay-rate differences between D® hys (t) — > 
X 5jL + 7r° and Dl hys (t) -> X 5 ,L + 7r° can be as large as 3%. Since 7^ s "°(l/r) ~ -7^ i7r °(l/r), 
it is possible to extract the rough size of Y' ~ — x D sin(20). Clearly the measurements of 
TZ^ sn \t), n^\t) and K* S7T (t), n^ Ln \t) are complementary to one another and can shed 
some light on the mixing parameters x D and y D as well as the possible new physics hidden in 
D° — D° mixing. 
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Note that the above discussions can be directly extended to neutral D decays to the final 
states like Kg,L + p°, Ks,l + a? and K SL + u, which occur through the same quark diagrams 
as D°/D° -> K s ,l + vr° (see fig. 2). Because anc i y(>) depend only upon the D° - D° 
and K° — K° mixing parameters, a sum over the above modes is possible, without any dilution 
effect on the signals of D° — D° mixing and CP violation, to increase the number of decay 
events in statistics. 

C. Final- state interactions in D — > KK 

Recently the CLEO Collaboration has searched for CP violation in neutral D decays to the 
CP eigenstates K + K~, Ks<p and Kstc . The confidence intervals (90%) on CP asymmetries 
in these three modes were found to be —0.020 < Akk < 0.180, —0.182 < Ak s </> < 0.126 and 
—0.067 < Ak s7T o < 0.031, respectively |10| . Although a definite signal of CP violation was not 



established from the data above, the possibility that these decays may accommodate CP asym- 
metries at the percent level could not be ruled out. In the following we shall concentrate on the 
final-state interactions in D°/D° — > K + K~ and K°K°, since they may affect the magnitudes 
of CP asymmetries significantly. The similar discussions can be extended to some other decay 
modes such as D°/D° — > -n + n~ and K s tt°. 

We begin with an isospin analysis of D° -> K + K~, D+ -> K + K° and D° -> K°K°. To 
do this, we assume that there is no mixture of D — > KK with other channels. In the language 
of quark diagrams [26|| , these modes can occur through both tree-level and penguin diagrams. 



However, such a naive description is problematic due to the presence of final-state rescattering 
effects 123- The final states KK may contain 1 = 1 and 1 = isospin configurations, and the 
overall decay amplitudes of D — > KK can be written as 

A + _ = (K + K-\H\D°) = \(A X + A ) , 

A 00 ee (K°K°\H\D ) = \ (Ai - Ao) , 

A +0 ee (K + K°\H\D + ) = A 1 , (4.31) 

where A\ and A are two isospin amplitudes. Clearly three decay amplitudes can form an 
isospin triangle in the complex plane: A + _ + A 00 = A +0 . Since the branching ratios of D — ► 
KK have been measured, one is able to determine the relevant isospin amplitudes from the 
relations above. For our purpose, we are more interested in the ratio of two isospin amplitudes: 
A /A 1 ee Z exp(i<£>). It is straightforward to obtain 

Z = y/2 + R 00 ) - 1 , cos^ = R +-- Ro ° , (4.32) 

where = \A^ /v4 + o| 2 and Rqo = |Aoo/^4+o| 2 are two observables. If the annihilation 

diagrams and penguin effects in D — » KK are negligible, then A^ , Aqq and A + q have a 
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common KM factor (i.e., V CS V* S ) from the dominant tree- level (spectator) quark transitions. In 
this case, <p is purely a strong phase shift and the magnitude of Z is independent of the KM 
matrix elements. 

Current experimental data give B(D° -> K+K~) = (4.54±0.29) x KT 3 , B(D° -> K°K°) = 
(1.1 ± 0.4) x 10" 3 and B(D + -> K + K°) = (7.8 ± 1.7) x 10" 3 0. The lifetimes of D° and D + 
mesons are r D0 = (0.415 ± 0.004) x 10~ 12 s and t d+ = (1.057 ± 0.015) x 10~ 12 s, respectively. In 
the neglect of small phase space differences of three decay modes, we obtain R + _ = 1.48 ±0.45 
and Rqq = 0.36 ± 0.22. The sizes of Z and tp can be solved from eq. (4.32) with the inputs 
of and Rqo, but there is large error propagation in this procedure, particularly for cos<^ 
which is bounded by unity. For simplicity and illustration, we plot the allowed regions of Z and 
cos (p in fig. 6. One can observe that 1.7 < Z < 2.0 and 0.3 < cosip < 1.0 (the central values 

of and i?oo lead to Z « 1.6 and cos<£> ~ 0.68). This implies that significant final-state 

interactions may exist in the processes D — > KK. 

The isospin amplitudes Ai and A can be expanded in terms of the tree-level and penguin 
transition amplitudes . Without loss of generality, we write f] 

Ai = A 1T exp [i(0 T + S 1T )) + v4i P exp[i(0 P + 5 1P )) , 

A = A 0T exp [i (0 T + S 0T )) + A 0P exp [i (0 P + 5 p)) , (4.33) 

where 0t and 0p are the overall weak phases of tree-level and penguin diagrams respectively, <5 n x 
and S Q p (with n=l,0) denote the corresponding strong phases. Hence Z and ip defined above are 
complicated functions of A n T, A n p, 0t, 0p, ^nT and 5 n j>. Since new physics may significantly 
affect the penguin amplitudes, direct CP violation is possible to appear in D — ► KK. A 
constraint on the 1 = 1 penguin contribution to D — > KK can be obtained by observing the 
decay-rate asymmetry between D + — > K + K° and D~ — > K~K°: 

\(K+K°\H\D + )\ 2 - \{K~ K°\H\D-)\ 2 
\(K+K°\H\D+)\ 2 + \(K-K°\H\D-)\ 2 

— 2AnAip sin(0p — 0t) sin(<5ip — 5n) 



(4.34) 



A\ T + A\ F + 2Ai T A 1P cos(0 P — T ) cos(5i P — 5 iT ) 

Note that the weak phase difference |0p — 0t| m &y be rather small within the standard model, 

but some sources of new physics (e.g., the existence of the fourth quark family or an iso-singlet 

up-type quark |2I|) can significantly enhance it through the breakdown of unitarity of the 3x3 

KM matrix in the penguin loops [0. The direct CP asymmetries in D°/D° —>■ K + K~ and 

K°K° contain both / = 1 and / = penguin contributions, and the latter can in principle be 

3 Here we have neglected the contributions of tree-level annihilation diagrams to D° — > K°K°, which involve 
both V CS V* S and V c dV* d . These two graph amplitudes are expected to have large cancellation with each other 
due to the GIM mechanism iBTtl. 
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distinguished from the former with the help of eq. (4.34). In practical experiments, Tr+k- 
and T k oko are cleanly detectable on the ip(3.77) resonance (see eq. (4.15) for illustration). 
If one wants to calculate the decay-rate asymmetries between — > K ± + K$,l or between 
D°/D° — > Ks t L + Ks,l, then the CP violation induced by K° — K° mixing in the final states 
has to be taken into account. 



It is also argued that inelastic final-state interactions may affect D — > KK ||27|| . This kind 
of effect is possible to yield a nonvanishing rate asymmetry between the charged D decays to 
K + K° and K~K°, even though the penguin contributions are negligibly small. To justify the 
role of penguin transitions and inelastic final-state interactions, one has to rely on the future 
data on direct CP asymmetries in the decay modes under discussion. 



5 Neutral D decays to non-CP eigenstates 

We proceed to consider the case that both D° and D° mesons decay to a common non-CP 
eigenstate. Most of such decay modes occur through quark transitions of the types c — > s{ud) 
and c — > d(us) or their flavor-conjugate counterparts, and the typical examples are the Cabibbo- 
allowed decay D° — > K~n + and the doubly Cabibbo-suppressed process D° — > K + n~. Because 
neutral D decays to K ±r K^ are of particular interest for the study of D° — D° mixing and DCSDs 
in charm physics, we shall concentrate on them in this section. Of course similar discussions 
can be extended to other non-CP eigenstates. 

Note that D° — > K ±r R T and their CP-conjugate processes take place only via the tree-level 
quark diagrams, on which no new physics can have significant effect [0, |2"3]| . Thus the four 
transition amplitudes are factorized as follows: 

A k -tt+ = (VcsV* d )T a exp(i5 a ) , A K+7T - = {V cd V* 8 )T b exp(i5 6 ) , 

A K +n- = (V* s V ud )T a exp(i5 a ) , A K - W + = (V* d V us )T b exp(i8 b ) , (5.1) 

where T a and T b denote the real (positive) hadronic matrix elements, 5 a and 5 b are the corre- 
sponding strong phases. Defining hxn = T b /T a and 5ktt = S b — S a , we obtain 

Px-n+ ~ Pk+tt- ~ -A 2 /iE- 7r exp(i^ 7r ) (5.2) 

to a good degree of accuracy, where A ~ 0.22 is the Wolfenstein parameter. In the factorization 
approximation, the magnitude of h^-K is expected to be of 0(1). The strong phase shift Sk-k 
vanishes only in the limit of SU(3) symmetry f3lfl . To fit the recent CLEO result for D° — > 
KV [|, which gives \p K -^+\ 2 = (0-77 ± 0.25 ± 0.25)%, one finds 5 K7T ~ 5° - 13° from a few 
phenomenological models |33|, j3"4| . Of course a larger value for 5k-k cannot be absolutely 
ruled out from current experimental data because of many uncertainties associated with the 
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empiric models used to analyze nonleptonic D decays. Finally, the expressions of Xk-w+ an d 
Xk+tt- are obtainable from eq. (5.2) as 



Xr-tt^ 

Xr+tt- 



-X 2 h 



Kn 



— X 2 hKir 



exp [i(5 K 7T + 20)] , 
exp [i{5 K ir ~ 20)] , 



(5.3) 



where we have used the notation of q/p given in eq. (2.2). 

A. Incoherent D decays to K ± tt t 

A lot of attention has been paid to the time distributions of incoherent D decays to K ± tt t 
(see, e.g., refs. |15|, ^3], |34|). In particular, Browder and Pakvasa have given a quite detailed 
analysis of the implications of CP violation and final-state interactions in the search for D° — 
D° mixing from D® hys (t) — > K + n~ and Dp hys (t) — > K~n + f34|. Our subsequent discussions 
are complementary to their work on three points: (a) the CP- violating asymmetry between 
fpw(f) — > K~it + and L>p hys (t) — ¥ is analyzed; (b) the different effects of x D and y D on 

D ph (i) — > K ± n T and D® hys (t) — > K ± tx t are explored in detail; and (c) the time-independent 
measurements of these decay modes are considered. 



Up to 0(x 2 D ), 0{y 2 D ) and 0(A 4 ) for every distinctive term, the decay rates of D to K ± it T 
can be directly read off from eqs. (2.22) and (2.23): 



R(D° phys (t) -»■ K-7T+) ex |^- w+ | 2 exp(-n)|l + AV, 



[VD COS i 5 Kn + 20) 



+x D sin (5 K7T + 20)] Tt - ~ (x 2 D - y 2 D ) T 2 t 2 } 



R(D° phyB (t) -> K+7T-) oc \A K - n+ \ 2 exp(-Tt) S^l + X 2 h K , 



[VD COS ( S Kn ~ 20) 



+x D sin (5 K7T - 20)] Tt - \ (x 2 D - y 2 D ) T 2 t 2 } 



(5.4) 



and 



R(D° phys (t) ^ K+7T-) oc \A K ^ + \ 2 exp(-Tt)^X 4 h 2 Kn + X 2 h K7 r 

-x D sin(^-20)]n + ^r 2 t 2 } , 

) oc \A k -kA 2 Q *v(-^t)^X 4 h 2 KlT + X 2 h K K 
-x D sin(^ + 20)]n + ^r 2 t 2 } , 



R(D° phys (t) ^ K- 



[VD COS ( 5 Krr - 20) 



[Vd cos ( S Ktt + 20) 



(5.5) 



where r and r have been presented in eq. (3.2). To probe CP violation and D° — D° mixing, 
the following two types of measurables can be analyzed in experiments: 
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(1) The CP- violating asymmetry 

W) ~ R{D hys{t) ^ K - n+) + R(D° phys (t) - X+O • 
Explicitly, we get 

A K Tz{t) ps —\ 2 h Kn A D cos(20) (y D cos 5^ + x D sin 

+ sin(20) (y D sin 5 Klx - x D cos 5 K7T )\ Tt , (5.7) 

where the observable has been defined before (see eq. (4.6)). One can see that Ak-k^) are 
composed of two sources of CP- violating effects, that in D° — D° mixing (proportional to Ap) 
and that from the interplay of decay and mixing (proportional to sin(2</>)). The magnitude of 
-A.Kn{t) is constrained by both the DCSD amplitude \ 2 hx-K and the D° — D mixing parameters 
x D and y D . Since A D < 1CT 2 is expected, Ax-nit) can reach the percent level only when the 
sin(20) term is significantly enhanced by new physics. In new physics scenarios with y D <C x D 
p3|, |34]|, we get Axnit) ~ ^^k-k x d sin(20) cos^^Ti as a safe approximation. Taking h KlT ps 1.8 
(implied by Ip^.^+I 2 ps 7.7 x 1CT 3 f9j), x D < 0.086 and 5kk > 5°, one finds the restriction 
A K7r (t) < 7.5 x 10- 3 n. 

(2) The combined decay rate 

R K *{t) = R(D° phys (t) - K+ir') + R(D° phys (t) -> K-n+) . (5.8) 

By use of eq. (5.5), we obtain 

R Kv (t) ex \A K - 7r+ \ 2 exp(-Tt){2\ 4 h 2 Kn + r D T 2 t 2 

+2X 2 hKn [cos(20) (y D cos 5k-k — x D sin 5k-k) 

— Ad sin(2(f)) (y D sin 5 + x D cos 5 k-k) T^} , (5.9) 

where r D = (r + f)/2 defined in eq. (3.3) has been used. The three terms of RKn{t), which have 
different time distributions, come respectively from DCSD, D° — D° mixing, and the interplay 
of these two effects. Thus the detection of Rxnit) can determine the D° — D° mixing rate 
r D and distinguish it from the DCSD contribution. If \<p\ is not large (e.g., in the standard 
model), the interference term will be dominated by cos(20)(?/ D cos 5^ — x D sin 8kk) due to the 
smallness of A^. In this case, the information about y D might be obtainable if the contribution 
of x D to the interference term is suppressed by small 8k-k- To justify the possible magnitude of 
0, however, one has to combine the measurements of Axnit) and R^if). 

Now we take a brief look at the time-independent decay rates of D phys /D phys — > K ± n T . 
With the help of eqs. (2.24) and (2.25), we obtain 



Ak-k — 



R(D% ys -» K-n+) - R(D° phys -> K + ix~) 
R(D° phys -> K-n+) + R(D° phys - 
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« -X 2 h 



K- 



A D cos(20) (y D cos 5 K7r + x D sin S K7l 



+ sin(20) (?/ D sin S KlT - x D cos S Kn )} 



(5.10) 



i.e., Aktt is approximately equal to the value of Ax-nit) at t — 1/T. Similarly, one can calculate 
another CP-violating asymmetry: 



Ak-k 



R{Dl 



phys 



K" 



7T 



R{D° 



phys 



K-7T- 



R(D° 



phys 



R(D° 



phys 



A 2 /i 



A D cos(20) (y D cos 5^ - x D sin 5^) 
+ sin(20) (?/ D sin 5 K n + x d cos5 k -k)] {A 4 /^ 
+r D + X 2 h KlT [cos(20) (y D cos S Kn - x D sin 5 Kw ) 
- A D sin(20) (y D sin 8 Kv + x D cos 5 Kn ) } 

Taking « and sin (20) « ±1 for example, we find that Ak-k m &y be significant: 

X 2 h K K (y D sin 5k* + x D cos S K7T ) 



(5.11) 



A 4 ^L + 



(5.12) 



D 



Due to the suppressed rates of -Dph ys — ¥ K + tt and D® hya ~* K 7r+ > however, the measurement 
of Ak-k will be a stiff experimental challenge. 

B. Coherent D decays to K ± tt t 

At the "0(3.77) and -0(4.14) resonances, the K ± n T events may come from the coherent 
decays of (D® h D^^c pahs. The flavor of one D meson decaying to K ± n T can be tagged 
by detecting the other D decaying to the semileptonic states I X^. The overall rates for such 
joint decay events, up to 0(xj- ) ), 0{y 2 D ) or 0(A 4 ) for every distinctive term, are obtainable from 
eqs. (2.26) and (2.27) as follows 4 : 



R(r,R-K + )- 
R(I + ,K + tt-)- 

R(r,K + n-)- 
R{l + ,K-7i + )- 



oc \AA 2 \A K - 



oc \A\\A 



K- 



oc |.Aj| 2 |.A K -- 7r -t 



(2 - x 2 D + y 2 D 
<^2-x 2 D + y 2 D 



2X 4 h 2 K7T +(x 2 D + y 2 D ) 
2X 4 h 2 Kn +(x 2 D + y 2 D 



(5.13) 



and 



R(r,K-n 



'■X 



\A l \ 2 \A K - 7r+ \ 2 {2-3{x 2 D -y 2 D ) 



*The formulas with the assumption of y D <C x D and \q/p\ = 1 have been given in rcf. [||. 
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[y D cos(5 K it + 20) + x D sin(<5x7r + 20)] \ , 



R{1 + ,K + ti-) + oc |^| 2 |A^+| 2 {2-3(. 



X 2 D - Vd) 



+ A\ 2 h Kn 



[Vd cos(5 K7r - 20) + x D sm(5 K7T - 20)] \ , 



i?(r,K + 7r-) + oc lAH^-^l 2 <2\*h 2 Kir + 3(a? D + y 2 D ) 



+ 4A 2 /i^ 



[j/d cos(5 X7r - 20) - x D sm(S K77 - 20)] \ , 



R(l + ,K-n+) + oc \A l \ 2 \A K - 7T+ \ 2 h\' i h 2 Kn + 3(x 2 D + y 2 D ) 



+ 4A 2 /i i ^ 7r 



cos(Sktt + 20) - x D sm(5 K 7T + 20)] 



(5.14) 



Some discussions about these results are in order. 



(1) To an excellent degree of accuracy, we have 

R(r,K-ir + )- « R(1 + ,K + ti-)_ . (5.15) 
The joint decay rates R(l T , K ± n T )_ can be normalized by R(l T , K T n^) _ , and the resultant 



rate difference or sum reads 



R(l-,K+n-)- 
R(1-,K-tt+)_ 



R(l + ,K-n+)_ 
R(l+,K+n-)_ 
R(l + ,K-n + )_ 



-2r D A D , 



2X 4 h 2 Kn + 2r 



D > 



(5.16) 



R(l-,K-n+)_ R(l+,K+n-)_ 
where r D and Ao have been given in eqs. (3.3) and (3.6), respectively. Observation of the 
CP- violating asymmetry S*! ^ may be practically impossible due to the smallness of and 
r D . However, is expected to be measurable at a r-charm factory running on the -0(3.77) 
resonance. As we shall show in the next subsection, r D can be extracted from the joint decay 
rates R(K + n~, K + tt~)- and R(K~n + , K + tt~)-, thus a comparison of this measurement with 
that for will separately determine the magnitudes of D° — D° mixing and DCSD. This 
idea is interesting on the point that the relevant measurements are time-independent and the 
involved decay modes are only D°/D° — > K ±/ k^ l . 



(2) It is easy to obtain the rate asymmetry 

R(l-,K-ir + )+ - R(l + ,K + n- 



2A 



R(l-,K-7i+) + + R(l+,K+7i-) + Kni (5,17) 

where Ar-k has been given in eq. (5.10). Normalizing the joint decay rates R(l T , K ± tt t ) + by 
R(P,K^n ± ) + , we get 

_ ( _) R(l~, K+n-) + R(l+,K-n+) + 



R(l-,K-n+). 



R(l+,K+ir- 
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« -6r D A D - A\ 2 h Kn A D cos(24>) (y D cos5 K7T - x D sm5 Kn ) 
-sin(20) (y D sm5 Kn + x D cos5 K x)} , 



s{+) ^ R(l-,K+n~) + + R(l+,K-7r 



R(l-,K-n+) + R(l+,K+7i-) + 
« 2A 4 /i^„. + 6r D + A\ 2 h Klv [cos(20) (y D cos5 K 7r - x D sm5 K 7r) 

-A D sm(2(p)(y D sinS Kn + x D cos5 Kw ) . (5.18) 

From eqs. (5.11), (5.16) and (5.18), one can see the following relation among and Ak-k- 

o(-) _ oo(-) 

*" K g^T^ • (5 - 19) 

This result could be tested if the data on all six measurables were available. 

(3) To give one a feeling of ballpark numbers to be expected, we roughly estimate the 
magnitudes of the above-mentioned observables by assuming A D = A D = and y D <C x D . 
Taking the semileptonic decay mode serving for flavor tagging to be D° —>■ K~e + u e , we have 
its branching ratio B(D° — > K~e + v e ) ~ 3.8% In addition, the current data give B(D° — > 
K-n + ) » 4.01% 0. Then R(l+,K+n-) ± and R(l-,K~n+)± are at the level 10~ 3 or so, while 
R{1~ i K + tc~)± and R(l + , K~tt + )± may be of the order 10~ 5 if we input x D ~ 0.06. Within the 
experimental capabilities of a r-charm factory, it is possible to measure the latter four decay 
rates to an acceptable degree of accuracy with about 10 7 D°D° events |TJj]. Furthermore, the 



upper bounds of the CP asymmetries Ak-k and can be obtained by use of the experimental 
results x D < 0.086 and Ip^-^+l 2 ~ 0.77% (I], [|. Taking cos^^ = 1 and sin(20) = ±1, 
we get \Aktt\ < 0.008 and \S+\ < 0.03. In the assumption of perfect detectors or 100% 
tagging efficiencies, one needs about 10 8 D°D° events to uncover \S ( + ] \ ~ 0.01 at the level 
of three standard deviations or to measure 21.4^1 ~ 0.005 in eq. (5.17) at the level of one 
standard deviation. Accumulation of so many events is of course a serious challenge to all types 
of experimental facilities for charm physics, but it should be achievable in the second-round 
experiments of a r-charm factory. 

C. Ratios of R(K ± 7r T , K^) c to R{K^, K^ir^c 

It has been pointed out that the coherent decays (L>p hys .Dp hys )c —* (K ± ti t )(K ± 7t t ) can be 
used to search for D° — D° mixing and to separate it from the DCSD effect f35fl . The relevant 
measurables are 

+ _ = R(K+n-,K+n-) c _+ = R(K-^,K~n+)c , . 

° ~ R(K-7T+,K+7T-)c ' ° ~ R(K-7T+,K+7T-) C " { ' 

Since in previous calculations the effects of CP violation or nonvanishing 5kh on were ne- 
glected, it is worth having a recalculation for these observables without special approximations. 
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By use of eqs. (2.26), (2.27) and (5.3), we obtain 

R(K-7r + ,K + 7r-)_ oc \A K - n+ \ 4 [2-x 2 D + y 2 D -AX 4 h 2 Kn cos(25 K7r ) 
R(K + n-,K-n+)_ oc \A K - n+ 1 4 [2 - x 2 D + y 2 D - AX 4 h 2 K7T cos(2^) 

oc \A K - w +\ 4 (x 2 D + y 2 D j 



and 



R(K-7i + ,K-n + 
R(K + ir-,K+n-)_ oc \A K - W+ 1 4 (x 2 D + y 2 D ) 



R(K-7i + ,K + n-) + oc \A K - n+ \ 4 {2-3x 2 D + 3y 2 D + 4\ 4 h 2 Kw cos(25 K7r ) 



+AX 2 h Kn 



AX 2 h 



Kw 



[y D cos (6 K „ + 2(j))+x D sin (5 KlT + 20)] 

[y D cos (5 X 7t - 20) + x D sin (5 Kw - 20)] 
R(K + n-,K-7i + ) + oc |A x -, + | 4 {2-34 + 3^ + 4A 4 ^ 7r cos(25 X7r ) 

cos (5xtt + 20) + x d sin (5 KlT + 20)] 



+4A 2 /i X7r 
+4X 2 h Kw 

R(K-7t + ,K-tv + ) + oc |A x -^+| 4 

+ 8\ 2 flKTT 



[y D cos (5 X 7t - 20) + z D sin (5 x?r - 20)] 

2 



3 (a& + y 2 D 



[y D cos (6 Kir + 20) - x D sin (5^ + 20)] 



R(K + 7r-,K + 7r-) + oc |A x -^ + | 4 3(x 2 D + ^ 



+ 8A 4 /i 



4 (.2 



Kit 



+ 8\ 2 h,Kir 



(5.21) 



cos (<5 K7r - 20) - x D sin (5 x?r - 20)] \ (5.22) 



up to O(a^), 0(y 2 D ) or 0(A 4 ). Clearly R(K-tt + , K + tt-) c ~ R(K + tt-, K-tv + ) c holds to an 
excellent degree of accuracy. As a consequence, the ratios r^ T are given by 



Xd + Vd 



x D + Vd 



(5.23) 



and 



3rl~ + AX 4 h 2 K7r + AX 2 h 



x [y D cos (5^ - 20) - x D sin (5^ - 20)] , 



-+ 



3rZ+ + 4X 4 h 2 Kn + 4X 2 h Kw 



x [y D cos (^tt + 20) - x D sin (5^ + 20)] 



(5.24) 
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One can see that and rZ + are approximately equivalent to r and f obtained in eq. (3.2). 
The difference between rt~ and rZ + measures CP violation in D° — D° mixing, and the sum 
of them amounts approximately to r D given in eq. (3.3). The DCSD effect on and rZ + is 
significant and non-negligible, but its magnitude can be isolated from the difference r^T — 3r+~ 
or rZ + — 3rZ + - In addition, we find 

r~ + - r+~ ^ 8X 2 h K7T A D cos(20) (y D cos 8 K * ~ %d sin $Kn) 

-sin(20) (y D sm5 K n + x D cos 5 K n)] ■ (5.25) 

Comparing this result with those derived in eq. (5.18), one gets 

rZ + - r+- « -2(6r D A D + S { + ] ) . (5.26) 



Such a CP- violating signal might be detectable at a r-charm factory running on the ■0(4.14) 
resonance. 

Although the above discussions concentrate only on D°/D° — > K ± 7r^ 1 , the similar results 
can be obtained for some other decay modes taking place via the same quark diagrams, such 
as D° — » K ± p T , K* ±r K T and their flavor-conjugate processes. All these channels are expected 
to have the same weak interactions, but their final-state interactions may be different from 
one another (e.g., 5 Kn ^ 5 Kp ). If the SU(3) breaking effects in D°/D° -> ( J ftT ± ,X* ± ) + 
(7r T , p T ,af, etc) are not so significant that all the strong phase shifts lie in the same quadrant 
as 8 Kir, then a sum over these modes is possible to increase the number of decay events in 
statistics, with few dilution effect on the signal of D° — D° mixing and CP violation. 



6 On CP- forbidden decays 

We now consider CP-forbidden transitions of the type 

(D° phys D° phys ) ± - (fJ 2 ) T , (6.1) 

where the D°D° pair with definite CP parity can be coherently produced on the ^(3.77) or 
■0(4.14) resonance, f\ and f 2 denote the CP eigenstates with the same or opposite CP parity. 
It is worth remarking that for such decay modes the CP-violating signals can be established 
by detecting the joint decay rates other than the decay-rate asymmetries. In practice, this 
implies that neither flavor-tagging for the initial D mesons nor time- dependent measurements 
of the whole decay processes are necessary. The joint decay rate P(/i, fzjc an d its analytical 
approximation have been presented in eqs. (2.21) and (2.26). For simplicity and illustration, 
here we concentrate mainly on the CP-forbidden decays (D phys D phys )- — > (/i/2)+, such as 
(A/ 2 ) + = (K+K-)(n+n-) and {K+ K~){K+ K~) . The case (P° hys P° hys ) + - (/i/ 2 )- will be 
briefly discussed by taking fi = Ksir and f 2 = Kltt for example. 



30 



By use of the quantities Uf, V/ and W/ denned in eq. (4.2), the joint decay rate R(fi, f 2 )- 
can be written as 

2 

x 



RU1J2)- oc I^J 2 !^! 2 (l + lA^I 2 ) (l + |A /2 | 2 ) 
1 , 1 



_ (i _ W h W f2 ) 
Ud 



X 



v 



(u fl u h + v fl v f2 ) 



(6.2) 



Here we assume f\ and f 2 to be two CP eigenstates with the same CP parity. CP conservation 
requires V h = V h = 0, U h = U h =0 and W fl = W/ 2 = ±1, then we get R(f u / 2 )_ = 0. Thus 
nonvanishing R(fi, / 2 )_ is a clean signal of CP violation. In the special case fi — f 2 = f, one 
finds 

«(/,/)- oc Mll + My^^-^M + V!) . (6.3) 

This result can be straightforwardly obtained from eq. (6.2) with the help of eq. (4.3). As 
discussed before, Uf is composed of the CP asymmetry in D° — D° mixing and that in the 
direct transition amplitudes of D decays, while V/ signifies the CP asymmetry induced by the 
interplay of decay and D° — D° mixing. Due to the smallness of Uf, V/, x D and y D , we believe 
that R(f, /)_ must be significantly suppressed. 

In comparison with R(f, /)_, the joint decay rate R(f, /) + is not CP-forbidden: 



R(f,f) 



+ 



(X 



'(l+IA/ry 

4y D 



p 


2 - 


q 





i + yh 



(l + Wff 



2 \2 



vh) 



Wf 



x 



D 



2 \2 



K 2 - V/ 2 ) 



Ax 



D 



2 \2 
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(6.4) 



Approximately, we obtain 



R(f,f)- „ (*2> + i/20(w? + y?) 



(6.5) 



R(f,f) + l + Wj-4y D W f ■ 

This relation can in principle be tested for / = K + K~ etc at the ■0(4.14) resonance in the 
second-round experiments of a r-charm factory, if the rate of D° — D° mixing is at the detectable 
level. 

In the neglect of CP violation in K° — K° mixing, the states K s n° and K L n° are two CP 
eigenstates with the opposite CP parity. Thus the process (Z)p hys Z)p hys ) + — * (K s tt°)(K l tt ) 
should be CP-forbidden. As a good approximation, we have |A^ i7r o| |Ax 57r o| and P Kl7T o ~ 
~Pk S tt° ( see ec f (4-22)). Then the joint decay rate with C = + turns out to be 
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x 



R(K L 7r°,K s n°) + oc \A Ks7T o\ 4 (l + \\ Ks7T o\ 2 ) 

i + yl 1 



x 



D 



;i - ylf 



(Uk sw o + V^ s7r o) 



(6.6) 
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Using the approximate results in eq. (4.24) and taking |e| ~ 0, we obtain a simpler expression 
for the equation above: 

R(K L 7i°,K s n°) + oc 6\A Ks7T o\ 4 (l + w)(x 2 D + y 2 D ) [A^ + sin 2 (20)] , (6.7) 

where w has been defined in eq. (3.3). In contrast, it is easy to check from eq. (6.3) that 

R(K L 7r°,K L n°)_ « R(K s ir°, K s n°)_ » \r{K l h\ K s ^) + . (6.8) 

Note that CP violation in D° — D° mixing (i.e., A D ) might be negligibly small, thus the 
dominant signal of CP violation in R(K s ir°, K s n )- or R(K L n°, K s n°) + could come from the 
mixing phase enhanced by new physics. In this sense, it is worthwhile to experimentally 
search for the above-mentioned CP-forbidden transitions. 



7 Summary 

To meet various delicate experiments in the near future at fixed target machines, P-meson fac- 
tories and r-charm factories, we have made a further study of the phenomenology of D° — D° 
mixing and CP violation in neutral P-meson decays. The generic formulas for the time- 
dependent and time-integrated decay rates of both coherent and incoherent D°D° events were 
derived, and their approximate expressions up to the accuracy of 0{x 2 D ) and 0{y 2 D ) were pre- 
sented. A variety of D° — D° mixing and CP- violating signals were analyzed in detail for 
neutral D decays to the semileptonic states, the nonleptonic CP eigenstates, the nonleptonic 
non-CP eigenstates and the CP-forbidden states. 

In particular, we have shown that it is possible to separately determine the magnitudes of 
x D and y D through precise measurements of the dilepton events of coherent D°D° decays on the 
'0(4.14) resonance at a r-charm factory. We gave a detailed analysis of D° — D° mixing signals 
and DCSD effects in the time-dependent and time-independent decays D°/D° — > K ± ir T . It is 
found that some constraints on x D and y D can be achieved in both fixed target and r-charm 
factory experiments, and the mixing and DCSD effects are distinguishable from each other. 
Taking CP violation and final-state interactions into account, we recalculated the joint decay 
rates of coherent D°D° pairs to (K ± -k t )(K ± 7t t ) , which are useful for the time-independent 
determination of r D and DCSD amplitudes. A special attention has been paid to the D° — D° 
mixing signals in the decay modes D°/D° — > Ks,l + 71-0 etc. We pointed out that a model- 
independent restriction on x D and y D should be obtainable from the time distributions of such 
decay modes. 

CP violation in D° — D° mixing can be well constrained in the semileptonic decays of 
coherent or incoherent D°D° events. In addition to this source of CP asymmetry, we have 
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shown that both the direct CP asymmetry in the transition amplitudes of D decays and the 
indirect CP asymmetry arising from the interplay of decay and D° — D° mixing can also 
manifest themselves in neutral D decays to hadronic CP eigenstates. These different CP- 
violating signals usually have different time distributions in the decay rates, thus they are 
possible to be distinguished from one another. In particular, direct CP violation can be cleanly 
probed in the coherent (P°P°)_ decays to a CP eigenstate plus a semileptonic state on the 
if) (3. 77) or ■0(4.14) resonance. For the decay modes with K° — K° mixing in the final states, 
however, the CP asymmetry induced by the mixing parameter e may be non- negligible and 
even dominant over the direct CP- violating signal from the charm quark transitions. Taking 
D — ► KK for example, we illustrated the significant effects of final-state interactions on CP 
violation. Different from those neutral D decays to CP eigenstates, D°/D° — ► K ± tt t are 
expected to have no direct CP asymmetries. Although the indirect CP-violating effects exist 
in such processes, they are suppressed to some extent by the DCSD amplitudes. We also 
discussed the CP-forbidden transitions on the 0(3.77) and 0(4.14) resonances. A search for 
the CP-forbidden modes like (PAD°)_ -> (tf+tf-) (tt+tT) and (D°D°) + -> (K s tt°)(K l tt ) is 
worthwhile in the future experiments of charm physics. 

Throughout our calculations CPT symmetry in the D° — D° mixing matrix has been 
assumed. Also the AQ = AC rule was assumed to hold in most cases, but the effects of 
AQ = —AC transitions on D° — D° mixing and CP violation were briefly discussed in sec- 
tion 3-B. Due to the smallness of x D and y D , it will be very difficult to accurately test the 
AQ = AC rule and CPT invariance in the D° — D° system. Recently Colladay and Kostelecky 
have studied a few possibilities to examine CPT symmetry in neutral D decays on the basis of 
the future fixed target and r-charm factory experiments |3"E|] . Considering this work and some 
other works on tests of CPT symmetry in the B° — B° system |37, 3"8| , we want to remark 



that one of the most sensitive signals for CPT violation or AQ = —AC transitions should be 
the nonvanishing asymmetry defined in eq. (3.5). However, one should keep in mind that 
Ajj might also come from the phase shifts of final-state electromagnetic interactions or the 
CP- violating contributions of non-standard electroweak models to the tree-level W- mediated 
semileptonic D decays. Another possible way to test CPT invariance in D° — D° mixing, which 
in principle works, is to measure the time distributions of opposite-sign dilepton events at an 
asymmetric r-charm factory (see Appendix B). 

Of course much more theoretical effort should be made to give reliable numerical predicitions 
for the magnitudes of various D° — D° mixing and CP- violating phenomena. 
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Appendix A 



This appendix is devoted to giving some generic formulas for the time-dependent D decays 
at an assumed asymmetric r-charm factory. The asymmetric e + e~ collisions just above the 
production threshold of (Dp hys /)p hys )c pairs will offer the possibility to measure the decay-time 
difference t_ = (i 2 — ti) between -Dp hys — > fi and -Dp hys — > f 2 . Usually it is difficult to measure 
the t + = (t 2 + ti) distribution in either linacs or storage rings, unless the bunch lengths are 
much shorter than the decay lengths [Q. Here we calculate the t_ distributions of joint decay 
rates starting from the master formula in eq. (2.19). For simplicity, we use t to denote t_ in 
the following. Integrating R(fi,ti, f2,t 2 )c over t + , we obtain the decay rates (for C = ±) as 



R(fi,f 2 ;t)- oc |A /l | 2 |A /2 | 2 ex P (-r|t|) x 

JM 2 + |C~| 2 ) cosh^H) - 2Re (CC-) sinh^H) 
^ )U 2 ~ IC-I 2 ) cos(x D n) + 21m (CC-) sm(x D Tt) 



(Al) 



and 



R(fi,f 2 ;t) + oc |A /l | 2 |A /2 | 2 exp(-r|t|) x 

'K + I 2 + IC 12 



ie + i 2 -ic + i 2 



x 2 D 



COSh(?/ D r|t| + (j)y) 

cos (x D r\t\ + 4> x ) 



2Re (C+Ch 



+ 



2Im(e;C+ 



x 2 D 



smh{y D T\t\ + 



sin(a; D r|t| + <f) x ) 



(A2) 



where the phase shifts cj) x and 4> y are defined by tan0 x = x D and tanh0 y = y D respectively. 
One can check that integrating R(fi, f 2 ;t) c over t, where t G (—00, +00), will lead to the 
time-independent decay rate R(fi, f 2 )c given in eq. (2.21). Eqs. (Al) and (A2) are two basic 
formulas for investigating coherent D°D° decays at asymmetric r-charm factories. 

Another possibility is to measure the time-integrated decay rates of (Dp hys Z)p hys )c with a 
proper time cut, which can sometimes increase the sizes of CP asymmetries [|L3[]. In practice, 
appropriate time cuts can also suppress background and improve statistic accuracy of signals. 
If the decay events in the time region t 6 [+t , +00) or t 6 (—00, — 1 ] are used, where t > 0, 
the respective decay rates can be defined by 



R(flj2,+t )c 

R(fi, f 2 ', —to)c 



-to 



R(fi,f 2 ;t) c dt, 



[ °R(fi,f 2 ;t)c dt. 



(A3) 
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By use of eqs. (Al) and (A2), we obtain 

£(/i,/ 2 ;±t )- oc \A h \ 2 \A h \ 2 exp(-Tt ) x 

i£-i 2 + lC- 12 



V 1 - ti> 
im 2 -ic-I 2 



Re (CC- 

cosh(y D n + y ) =F — i ' smh(y D Tt Q + Yy) 

i-vb 



Im ( £L C- 

cos (x^Tto + E ) ± — 7 sin(x £) rt + 

1 + x 2 D 



and 



R(f 1 ,f 2 ;±t )+ oc |%| 2 |A /2 | 2 exp(-n ) x 

^+I 2 + IC 12 



2(1- 2/1,) 

ie+i 2 -ic + i 2 

2(1 + 4) 



Re f£* (T ) 

cosh(?/ £) rto + 20 y ) — ±_!— sinh(y D n + 2<j) y ) 

1 ~ Vd 
Im(aC+) 

cos {x D rt + 2(p x ) H — — 2 — sm(x D rt + 2(f> x ) 

1 + X D 



It is easy to check that 



(A4) 



(A5) 



£(/i,/ a ;+o) c + i2(/i,/ 2 ;-o)c = R{hJ 2 )c. 



(A6) 



One can observe that in / 2 ; ±io)c different terms are sensitive to the time cut to hi 

different ways. Thus it is possible to enhance a CP-violating term (and suppress the others) 
via a suitable cut to. 



Appendix B 

In this appendix we take a brief look at the possible effect of CPT violation in D° — D° 
mixing on the decay rates of semileptonic D decays. For simplicity, we assume the AQ = AC 
rule and direct CPT invariance in D decays to hold exactly. We also assume the absence of 
final-state electromagnetic interactions and other sources of new physics that could affect the 
tree- level VT-mediated D decays. Due to the presence of CPT violation, the mass eigenstates 
\D L ) and \Du) can now be expressed as 



\D L ) = cos^ p\D°) + sin^ q\D°) 



9 9 
\D n ) = sin- p\D°) - cos- q\D } 



;bi) 



where 9 is in general complex. Note that CPT invariance requires cos 9 = 0, while CP conser- 
vation requires both cos# = and p = q = 1 j|(J. Taking 9 = 7r/2, i.e., CPT symmetry, one 



can reproduce eq. (2.1) from eq. (Bl). The proper-time evolution of an initially (t = 0) pure 
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D° or D° turns out to be 



\D° phys (t)) = [g + (t)+9-(t) cos 9} \D°) + l[g_(t)sm8}\D°) , 
\D° phys (t)) = l9 + (t) ~ 9-(t) cose] \D°) + V - [g_(t) sin0] \D°) , 



(B2) 



where g±(t) have been given in eq. (2.6). 

Starting from eq. (B2), one can calculate the CP asymmetry A D defined in eq. (3.5) for 
semileptonic D transitions. We find 

-r 2x D a:Im(cos#) + 2?/ D Re(cos#) . . 

D = (1 + a) + (l-a)|cosfl| 2 ' ( } 

where a = (1 — y 2 D )/(l + x 2 D ) has been defined before. Clearly = 0, if there is no CPT 
violation in D° — D° mixing (i.e., cos# = 0). Since | cos#| must be a small quantity, the | cos#| 2 
term in the denominatior of is negligible. Anyway observation of the signal Ap will be 
greatly difficult in practice, since its magnitude is suppressed by both the small mixing rate 
and the small CPT asymmetry. 

Next let us assume the experimental scenario to be an asymmetric r-charm factory, in 
which D°D° pairs can be coherently produced at the ip(3.77) or ■0(4.14) resonance and the 
time-dependent measurements of their decays are available. To probe possible CPT violation 
in D° — D° mixing, we consider the case that one D meson decays to the semileptonic state 
e ± Xf at (proper) time t e and the other to the semileptonic state fJ^X^ at t^. The joint decay 
rate for having such an event can be given as a function of the decay-time difference t = t^ — t e . 
For simplicity and definition, we choose t > by convention. This implies that e ± Xf events 
may serve for flavor-tagging of fJ^X^ events. After a lengthy calculation, we obtain 



R(e ± X^,^X^;t)- oc |A e | 2 |A M | 2 exp(-n) [cosh(y D n) + cos(x D Tt) 

± 2Re(cos#)sinh(y D rt) ± 2Im(cos 6) s\n(x D Yt)} , 



(B4) 



and 



R(<*X?,^X};t) + oc \A e \M»? exp(-rt) j + ^ + + ^ 

2 1 cos#| 

± [cos(6 + ojJ) exp(+y D n) - cos(0 + u + + x D Yt)\ (B5) 

/4 + (2-y D ) 2 

21 cos #1 1 
=F , = [cos(0 — uj + ) exp(—y D Tt) — cos(0 — u + — x D Tt)] > , 

Jx* D + (2 + y D f J 

where <p x and <p y have been defined in Appendix A, and the phase shifts uj± and are defined 
by tancj ± = x D /{2 ±?/£>) and tan© = Im(cos#)/Re(cos#) respectively. In obtaining eqs. (B4) 



36 



and (B5), we have neglected those higher-order terms of cos 8. It is clear that the opposite-sign 
dilepton events i?(/ ± X[ F , PX t ; t)c cannot be used to explore possible CPT violation in D°—D° 
mixing, because the time order of l + and l~ is hardly distinguishable in practical experiments. 
In addition, the signal of CPT violation cannot manifest itself in the time- integrated decay 
rates of (-Dphys-^phys)c — ► ( e± ^ =F )(A i=F ^ ± )) as obviously shown by the equations above. That is 
why we need an asymmetric r-charm factory to test CPT symmetry in D° — D° mixing. 

Of course, CPT violation can appear in many other decay modes of neutral D mesons. The 
semileptonic processes discussed above are more attractive to us for the study of CPT violation, 
since they do not involve CP asymmetry in D° — D° mixing (measured by \q/p\ ^ 1) and other 
CP-violating signals. In general, however, both direct and indirect CPT asymmetries as well 
as AQ = —AC transitions (and other sources of new physics) are possible to affect the decay 
modes in question ET . 



37 



References 

[1] Particle Data Group, L. Montanet et al, Phys. Rev. D 50, 1173 (1994). 

[2] E691 Collaboration, J.C. Anjos et al, Phys. Rev. Lett. 60, 1239 (1988). 

[3] H.Y. Cheng, Phys. Rev. D 26, 143 (1982); A. Datta and D. Kumbhakar, Z. Phys. C 27, 
515 (1985). 

[4] L. Wolfenstein, Phys. Lett. B 164, 170 (1985); T.A. Kaeding, Phys. Lett. B 357, 151 
(1995). 

[5] J.F. Donoghue, E. Golowich, B.R. Holstein, and J. Trampetic, Phys. Rev. D 33, 179 
(1986); G. Burdman, Report No. Fermilab/Conf-95/281-T (1995). 

[6] H. Georgi, Phys. Lett. B 297, 353 (1992); T. Ohl, G. Ricciardi, and E.H. Simmons, Nucl. 
Phys. B 403, 605 (1993). 

[7] LI. Bigi and A.I. Sanda, Phys. Lett. B 171, 320 (1986). 

[8] I.I. Bigi, in the Proceedings of the Tau-Charm Factory Workshop, SLAC, May 1989 (edited 
by L.V. Beers), p. 169. 

[9] CLEO Collaboration, D. Cinabro et al, Phys. Rev. Lett. 72, 406 (1994). 

[10] CLEO Collaboration, J. Bartelt et al, Phys. Rev. D 52, 4860 (1995); E687 Collaboration, 
PL. Frabetti et al, Phys. Rev. D 50, R2953 (1994). 

[11] D.M. Kaplan, in the Proceedings of the Charm 2000 Workshop, Fermilab, June 1994 
(edited by D.M. Kaplan and S. Kwan), p. 35; Report No. IIT-HEP-96/2 (1996). 

[12] J.R. Fry and T. Ruf, in the Proceedings of the Third Workshop on the Tau-Charm Factory, 
Marbella, June 1993 (edited by J. Kirkby), p. 387. 

[13] Z.Z. Xing, Phys. Rev. D 53, 204 (1996). 

[14] Z.Z. Xing, Phys. Lett. B 372, 317 (1996). 

[15] T. Liu, Report No. Princeton/HEP/95-6 (1995). 

[16] Z.Z. Xing, Phys. Lett. B 379, 257 (1996). 

[17] G. Gladding, in the Proceedings of the Tau-Charm Factory Workshop, SLAC, May 1989 
(edited by L.V. Beers), p. 169; U. Karshon, in the same proceedings, p. 706. 

[18] D. Du, Phys. Rev. D 34, 3428 (1986). 

38 



[19] I thank A.I. Sanda for useful discussions about this aspect. 

[20] See, e.g., LI. Bigi, F. Gabbiani, and A. Masiero, Z. Phys. C 48, 633 (1990). 

[21] Z.Z. Xing, Phys. Lett. B 371, 310 (1996). 

[22] J.L. Hewett, in the Proceedings of the Tau-Charm Factory in the Era of S-Factories and 
CESR, SLAC, August 1994 (edited by L.V. Beers and M.L. Perl), p.5. 

[23] G. Blaylock, A. Seiden, and Y. Nir, Phys. Lett. B 355, 555 (1995). 

[24] I.I. Bigi and H. Yamamoto, Phys. Lett. B 349, 363 (1995). 

[25] Z.Z. Xing, Phys. Lett. B 353, 313 (1995). 

[26] Z.Z. Xing, Nuovo Cimento A 108, 1143 (1995). 

[27] XY. Pham, Phys. Lett. B 193, 331 (1987); A.N. Kamal, Int. J. Mod. Phys. A 7, 3515 
(1992); A.N. Kamal and T.N. Pham, Phys. Rev. D 50, R1832 (1994). 

[28] W.F. Palmer and Y.L. Wu, Phys. Lett. B 350, 245 (1995). 

[29] Y. Nir, Report No. WIS-95-28-PH (unpublished); G.C. Branco, PA. Parada, and M.N. 
Rebelo, Phys. Rev. D 52, 4217 (1995). 

[30] Z.Z. Xing, Nucl. Phys. B (Proc. Suppl.) 50, 24 (1996). 

[31] L. Wolfenstein, Phys. Rev. Lett. 75, 2460 (1995). 

[32] L.L. Chau and HY. Cheng, Phys. Lett. B 333, 514 (1994). 

[33] F. Buccella et al, Phys. Rev. D 51, 3478 (1995). 

[34] T.E. Browder and S. Pakvasa, Report No. UH 511-828-95 (1995). 

[35] H. Yamamoto, Report No. CALT-68-1318, 1985 (unpublished). 

[36] D. Colladay and V.A. Kostelecky, Phys. Rev. D 52, 6224 (1995). 

[37] M. Kobayashi and A.I. Sanda, Phys. Rev. Lett. 69, 3139 (1992). 

[38] Z.Z. Xing, Phys. Rev. D 50, R2957 (1994). 

[39] G.J. Feldman et al, in the Proceedings of High Energy Physics in the 1990s, Snowmass, 
Colorado, 1988, edited by S. Jensen (World Scientific, Singapore, 1988), p. 561; K. Berkel- 
man et al, Report No. CLNS 91-1050, 1991 (unpublished). 

[40] T.D. Lee and C.S. Wu, Annu. Rev. Nucl. Sci. 16, 471 (1966). 

[41] Z.Z. Xing, unpublished note (1995). 



39 



0-8 I 1 1 1 1 1 1 1 r 




-0.8 1 1 1 1 1 1 1 1 1 1 

0.02 0.04 0.06 0.08 

x D 

Figure 1: The illustrative plot for the changes of 7 with x D , where the restriction 10~ 4 < r < 
3.7 x 10~ 3 has been used. 





Figure 2: Quark diagrams for the Cabibbo-allowed decay D° — > K 0r K° and the doubly Cabibbo- 
suppressed decay D° — > K°ir . 
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Figure 3: The illustrative plot for the time distributions of S + S7r (t) and S + L7r (t), where 
y D ~ 0.08 and (f> m have been taken. 
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Figure 4: The illustrative plot for the time distribution of Aisit), where y D « 0.08 and 
| cos(20)| Pi 1 have been taken. 
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Figure 5: The illustrative plot for the time distributions of s5 s7r (t) and s5 L7r (i), where 
x D Pi 0.08, ~ 7r/4 and Pi have been taken. 
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Figure 6: Possible regions of Z and cos 99 for I? — > KK (i.e., the quadrangle (abed)) allowed 
by current data: R + _ = 1.48 ± 0.45 and R m = 0.36 ± 0.22. 
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